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Abstract 

In this article, we construct a representation formula for stochastic B-series evaluated in a B-series. This formula is 
used to give for the first time the order conditions of implicit Taylor methods in terms of rooted trees. Finally, as an 
example we apply these order conditions to derive in a simple manner a family of strong order 1 .5 Taylor methods 
applicable to Ito SDEs. 
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1. Introduction 

Taylor methods have for a long time been a common choice for solving stochastic differential equations (SDEs), 
and with increased use of automatic differentiation techniques, their popularity will hardly subside. The Taylor ex- 
pansions of the exact solutions of SDEs, from which the Taylor methods are derived, can take one out of two forms: 
either as Wagner-Platen series [10, 11] or as B-series [6, 1, 2, 13, 15], the relation between the two series has been 
demonstrated in [8]. In the present paper, we focus on B-series. In short, the exact solution X{t) evaluated after one 
step starting at (fo, xq) can be expressed in terms of a B-series: 



X(to + h) = B{^p, XQ- h)^Yj ^(^) ■ ^(^)('') ■ ^(^)(-^o), 



where a is a combinatorial term, (/; is a stochastic integral, and F is some differential. These series will be thoroughly 
explained in Subsection 2.1. Implicit Taylor methods have been derived e.g. in [11, 16]. Written in terms of B-series, 
they take the form 

i-i = B(<1),,, xo; h) + B(0,„„ Yv,h). (1) 
To illustrate this, consider the semi-implicit Milstein scheme applied to the Ito SDE 

dX = f{X)dt + g{X)dW 

which is given by 

Fi = xo -h hf{Yi) + /(i)g(xo) + I(i,i)(g'g){xa). 
The Milstein method takes the form ( 1 ) with 

B(<1>„, xo; h)^ XQ + /(i)g(xo) -i- /(i,i)(g'g)(xo) 

and 

B{^i,„,Yuh)^hf(Yi). 
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To compare with the exact solution, the numerical solution Fi of (1) has to be expressed in terms of B-series as well. 
Here, this is done in two steps. In Subsection 2.2 we prove a quite general result, well known from deterministic 
theory: A B-series evaluated in a B-series is again a B-series, written as 

B{(t)y, B{<p„ xo; h); h) = B(,p, o <Py, xq; h). (2) 

This result is not only useful for implicit Taylor methods, but also for composition- and splitting methods. The rule 
can be used with both Stratonovich and Ito calculus. In Section 3, we present the B-series for the implicit Taylor 
methods, and thereby their order conditions. Finally, in Section 4, a family of implicit Taylor methods is developed 
and the results are confirmed by numerical experiments. 



2. Some notation, definitions and preUminary results 

In this section we introduce some necessary notation and provide a few definitions and prehminary results that 
will be used later. Let (D, f) be a probabihty space. We denote by {X{ty}tei the stochastic process which is the 
solution of a <i-dimensional SDE defined by 

m 

dX{t) = go(X(t))dt + 2 glim) * dWiit), X(to) = xo, (3) 
1=1 

with an m-dimensional Wiener process (W(f)),>o and / = [to, T]. As usual, (3) is construed as abbreviation of 

Xit) = xo+ goiXis))ds + Y, gi(Xis))-kdW,is). (4) 

The integral w. r. t. the Wiener process has to be interpreted e. g. as Ito integral with -kdWiis) = dWiis) or as 
Stratonovich integral with -kdWiis) = odWi(s). We assume that the Borel-measurable coefficients gi : M.'^ M.^ 
are sufficiently differentiable and that the SDE (4) has a unique solution. 

To simphfy the presentation, we define Wois) = s, so that (4) can be written as 

m „t 

X{t) = xo + Yj gi(X{s))*dWi{s). (5) 

1=0 ^'o 



In the following we denote by H a set of families of measurable mappings, 

H := [{(p{h)}h>o : <fi{h) : Q ^ M is ^ - S-measurable V/j > 0}. 



2.1. Stochastic B-series 

B-series for deterministic ODEs were introduced by Butcher [5]. Today such series appear as a fundamental 
tool to do local error analysis on a wide range of problems. B-series for SDEs and their numerical solution by 
stochastic Runge-Kutta methods have been developed by Burrage and Burrage [ 1 , 3] to study strong convergence in the 
Stratonovich case, by Komori, Mitsui and Sugiura [ 1 3] and Komori [ 1 2] to study weak convergence in the Stratonovich 
case and by RoBler [14, 15] to study weak convergence in both the Ito and the Stratonovich case. However, the 
distinction between the Ito and the Stratonovich integrals only depends on the definition of the integrals, not on 
how the B-series are constructed. Similarly, the distinction between weak and strong convergence only depends on 
the definition of the local error. Thus, we presented in [6] a uniform and self-contained theory for the construction 
of stochastic B-series for the exact solution of SDEs and its numerical approximation by stochastic Runge-Kutta 
methods. In this subsection, we review results given there for the solution of (5). In Subsection 2.2 we prove the 
composition rule (2). 

Following the idea of Burrage and Burrage we introduce the set of colored, rooted trees related to the SDE (4), as 
well as the elementary differentials associated with each of these trees. 
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Defuiition 1 (Trees). The set ofm + \-colored, rooted trees 

r = {0} u To u Ti u • • ■ u r„ 

is recursively defined as follows: 

a) The graph = [0]/ with only one vertex of color I belongs to Ti. 

Let T - [ti,t2, . . .,T^]i be the tree formed by joining the subtrees ti,T2,...,Tk each by a single branch to a common 

root of color I. 

b) //ti,T2,...,t^ e T them = [n, T2, . . . , T;,]; e T/. 

Thus, Ti is the set of trees with an /-colored root, and T is the union of these sets. As usual, the trees are not 
ordered, i. e. the order of subtrees is not significant. 

Definition 2 (Elementary differentials). For a tree t & T the elementary differential is a mapping F(j) : M"' — » M'' 
defined recursively by 

a) F(0)(xo) = xo, 

b) F{*i){xo) = giixo), 

c) IfT = [n, T2, . . . , rj/ e Ti then 

F{T)iXo) = gf\xo)iF{T,){Xo), F(T2)(Xo), . . . , F(Tj(Xo)). 

As will be shown in the following, both the exact and the numerical solutions can formally be written in terms of 
B-series. 

Definition 3 (B-series). A (stochastic) B-series is a formal series of the form 

B{<f>, xo; h) = Yj «W ■ ^We^) • F(t)(xo), 

where : T — > H and a : T Qis given by 

1 " 

a(0) = l, «(•/) = !, a(T=[Ti,...,TJO=-i-n rfl^^'^^ 

ri\r2\ ■ ■ ■ rg\ i.J 

where r\,r2, ■ ■ ■ ,rq count equal trees among t\,T2, ■ ■ ■ ,Tk. 

Note that a(T) is the inverse of the order of the automorphism group of t. To simphfy the presentation, in the 
following we assume that all elementary differentials exist and all considered B-series converge. Otherwise, one has 
to consider truncated B-series and discuss the remainder term. 

The next lemma proves that if Y{h) can be written as a B-series, then f(Y(h)) can be written as a similar series, 
where the sum is taken over trees with a root of color / and subtrees in T. The lemma is fundamental for deriving 
B-series for the exact and the numerical solution. It can also be used for deriving weak convergence results. 

Lemma 1. If Y(h) = B{^, xq; h) with 0(0) = 1 is some B-series and f e C°°(R'', M^) then f{Y{h)) can be written as a 
formal series of the form 

f{Ym - Yj • 'I'Mih) ■ G(m)(xo), (6) 
where Uf is a set of trees derived from T by 
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a) [0]/ e Uf, andifT\,T2,...,TK e T then [ti,T2,...,tJ/ e Uf, 



b) cm /)(xo) = f{xo) and G(m = [n , . . . , rj /)(xo) = /W(xo)(F(ti)(xo), . . . , F(tJ(xo)), 



1 * 

c) my) = 1 ^i^'A" = [Ti, . . . ,rj/) = ^ , fj a(T,), 

where ri,r2, . . . ,rq count equal trees among ti, ti, . . . , Tk, 

d) <i'^([0]/) = 1 an^; ^ir^(M = [Ti, . . .,T,-\f)Qi) = \Yj=x <l>(jj)Qi). 
Proof. Writing Y{h) as a B-series, we have 



CO ^ / 

°° 1 ' f " 



«t=l {ri,T2,...Tj€r\{0) 



/''^(xo)(F(ti)(xo),...,F(tJ(xo)) 



where the last sum is taken over all possible unordered combinations of k trees in T . For each set of trees ti , T2, . . . , T,f g 
T we assign a m = [ti, T2, . . . , t^]/ e f//. The theorem is now proved by comparing term by term with (6). □ 

When Lemma 1 is applied to the functions gi on the right hand side of (5) we get the following result: If Y{h) = 
B{(p, xq; h) then 

gmh)) = 2 a(T) • 0;(T)(/j) ■ F(r)(xo) (7) 



T€T, 



in which 



1 ifT=.,, 

K 



Theorem 2. The solution X{to + h) of (5) can be written as a B-series B{<p, xq; h) with 



m = h 



(f(T = [ti,...,tJ/)(/i) 



Cp(Tj)(s)'kdWl(s). 



Proof. Write the exact solution as some B-series X(to + h) = B((p, xq; h) with ^(0) = 1. By (7) the SDE (5) can be 
written as 

2 a(T) • ipirm ■ F{T)(xo) = + J] J] a(T) • ^;(r)(i) ★ dWi{s) ■ F(t)(xo). 

reT 1=0 reT, ^0 

Comparing term by term we get 

(p(T)ih)= f <p'(T)is)-kdWi(s) for r e T,, / = 0, 1, . . 
Jo 

The proof is completed by induction on r. □ 

The same result is given for the Stratonovich case in [3, 12J, but it clearly also applies to the Ito case. 

The definition of the order of the tree, p(t), is motivated by the fact that E Wi{h)^ = ft for / > 1 and Wo(ft) = h: 



Definition 4 (Order). The order of a tree r eT respectively u e Uf is defined by 

K 

p(0) = 0, p(M = [Ti , . . . , Tj f) = Y, P(rd, 

and 



i=\ 



P(t = [n, . . . , rJz) = J^p(Ti) + I 1 



(=1 



1 forl = 0, 
otherwise. 



We conclude this subsection with an example. 
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Example 



we get 



P(t) = 3, a(T) = -, 



f'(T)(xo)=g;,'tei,g'i'te2,^2))(xo) 



and 



<p{T)ih) = £ WiiSi) 1^ ' W2(S2f * dWi(S2)ySi 

4/(2,2,1,1,0) + 2/(2,1,2,1,0) + 2/(1,2,2,1,0) (Stratonovick), 

4/(2,2,1,1,0) + 2/(2,1,2,1,0) + 2/(1,2,2,1,0) 
+2/(0,1,1,0) + 2/(2,2,0,0) + /(1,0,1,0) + /(0,0,0) (ltd). 

2.2. Composition ofB-series 

In this section, we will show that, similar to the deterministic case, B-series evaluated in a B-series can be written 
as a B-series, or more specific: 

Bicf,y, B(4>„ xo; h); h) = B{^, o cPy, xq; h), (8) 

and we will define the composition operator <l>y 

For this purpose, we will need the decomposition of a tree t e T into a subtree d- sharing the root with t, and a 
remainder multiset o) of trees left over when & is removed from t. This forms a triple (r, co). We also include the 
empty tree as a possible subtree, in which case the triple becomes (r, 0, {t}). 

Example 2. Two examples of such triples are 




'W /y.i.,.,.,! 



corresponding to 



and 




or equivalently 



corresponding to 



T = 



5 



From this example we observe that the same pair of a tree r and a subtree & may have more than one remainder 
muhiset. We then define S T(t) as the set of all possible subtrees of t together with the corresponding remainder 
multiset ca, that is, for each t e r\0 we have 

5 r(.,) = {(0, {•;}), (•/,{m 

5r(T= [ti,...,tJ/) = : = [j?i,...,!?J/, w = {aJi,...,aj^}, 

{■&i,COd€ST{Ti), J= 1,...,^}U(0,{t}). 

Further, for each § € T define 

A(i?) = {(t,w) : e5T(T)}. 

At last, we will use that 

{(T,i^,aj) : &€T,iT,oj)€Ai^)} = [iT,^,oj) : r e T, w) e 5 r(T)} . (9) 
And with these preparations, we are ready to prove (8). From Definition 3 we have 

Bi<Py, BiiP,, xo; h); h) = Yj '^^^^ ' ^/^)(^) ' FimBi<p„ xq; h)). (10) 

We are now looking for an expression for F(i?) of the form 

FmB(<p,,xo;h))= Yj 7(T,&,oj)^-UT,oj)(h)-F(T)(xo). (11) 

In this case, we insert (11) into (10) and apply (9) to obtain (8) with 

{<Px°<Py)(T)^ Yj 7(T,&,aj)<pymMr,(o). (12) 

To find recursive formulas for the still unknown quantities y and we start with the trivial case i? = 0. We have 
A(0) = {(t, {t}) : reriand 

F(0)(B,) = B, 

where, for simphcity, we use the abbreviation B^ = B(<px, xo;h). Thus we can choose 7(t, 0, {t}) = 1 and ^^(t, {t}) - 
</>xiT)- For 1^ = we have A(«;) = {(r, co) : t = [ti, . . . , r^];, (o = {ti, . . . , r^} }. By Lenmia 1 and Definitions 1, 2, 
and 3 we get 

K 

F{»i){Bx) = giiBx) = J] aiT) ■ Y] (pxirdih) ■ F(t)(xo), 

reT, (=1 

SO we can choose y(T, &,a>) = I and 4>x(t, cS) - R/Li 4>x(Ti). For the more general case, let = . . . , so that 

FimBd = 8f"\Bx) {Fm(Bx), F(§,J(Bx)) . 

Apply Lenmia 1 on gf*\Bx), use the induction hypothesis (11) on each F{'&i){Bx) and the multUinearity of gf*^ to 
obtain 



(13) 



^|'^°^'^*>(xo) {F{6x){xo\ F(6J(xo), F(ti)(xo), . . . , F(t,„)(xo)) , 



where each = [Si,. . . ,5,fJ^«#. The last term can be recognized as the elementary diff'erential F(t)(xo) for r 
[^1, . . . , Ti, . . . , T;,^]/. Choosing 



^x(t, aS) = 



V !=1 / ;=1 
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we obtain by induction 



Seoj 



To complete the argument we have to find the expression for y(T, w) for this (t, i?, oj) triple. For this purpose, 
write?? = [d-\\...,d-''g']i, in the sense that ■& is composed by q different trees, each appearing r, times. Further, let 
Tik £ {ti,.. .,Tk] be the distinct trees with multiphcity r,jt, k = I, pi, of which the &i are subtrees, ^^li = ''(• 
Finally, let 5, e be the distinct trees with multiplicity Si,i= I,..., po, of the remainder multiset which are directly 
connected to the root of r, 2^=i = ko- Then, t can be written as 



■ ' ^qpq I' 



(14) 



where the rightmost expression above indicates that t is composed by Q different trees each appearing Ri times. Now 
the terms corresponding to each of the subtrees [rj^j', . . .,TjJ']; appear exactly ri\/{ra \ ■ • -fip,!) times in the sum (13). 
Using Picoo) from Lennma 1 the combinatorial term in (12) becomes 



7(t, oj) 



;n''<«jin;^)in^<-'.»"".)5i* 

• i=i )\i=i '^'i- '^'P'-j\i=i "^^f 



aid-) [si\---Spg 
From the definition of a for r and & we obtain 

y(T, 1?, cj) = 



Rll -Rnl 



Let us illustrate the use of this formula by an example: 
Examples. LetT= ■ Then for 

■& = we have u) = {•j, •;■} and y{j, §,(o) = I, 



(15) 



and for 

Next consider the tree 

One possible partition of this tree is 
Then (14) becomes 



we have oj = {•j] and y(T, ■&,o)) = 2. 



(16) 



T = 




(17) 



J' J' *] 7' *i 



and we get kq = 2, po = 1, si = 2, % = 2, rn = ri2 = I, Q = 2 andRi = 3, R2 = 1, so for the partition (17) ofr we 
get 

7(t, &,oj) = 3. 



For the partition 



(18) 



7 



(14) becomes 



7' *] i' *] 



giving K(j = 3, pq - \ so si — 3. Further — 1, and Q, Ri andR2 are as above. For the subtree we have to 

use the partition (16), thus for the partition (18) we get 

7(t, §, aS) = 2. 

We can now state the main result of this section: 

Theorem 3 (Composition of B-series). Let (px^fpy '■ T E and (pxi^) = 1- Then the B-series B{(j)y, •; h) evaluated 
at B((f>jc, xo; h) is again a B-series, 

B{4>y, Bi<px, xo; h); h) = Bi^^ ° <Py, xq; h), 

where 

{<Px ° <Py) (T) = Yj ^' '^^ '^y^^^ n ^^^^^ 

(i>,w)€Sr(T) V Seu 

with 7(0, 0, {0}) = 1 and yir, a>) given by (15) for Ti= 0. 

Remark 1. The composition rule cpx ° (/>y can also be expressed in terms of ordered trees, as in the deterministic case, 
see e. g. [9]. The combinatorial term y{T, oS) is then the number of equal terms appearing in the sum over ordered 
subtrees. 

Example 4. For example 

{<l>x ° <Py) ('M^'j = <Pym<Px ("^'^^''j + 'l>y(»j)M»l)^M»j) 
+ 24>y (/>xi»l)(/>xi»j) + (Py (I -j 4>x{*lf 
+ 203- V/ h*x(»/) + 0y V/ h*x(»j) + 0). 



and 



3. Application to implicit Taylor methods 

We will now see how the composition rule will give the order conditions for implicit Taylor methods. One step of 
this method is given by 

Fi = Bi(^ex, xo; h) + BiOi„, Fi; h) (19) 

with <l)e.j(0) = 1 and <I)„„(0) = 0. If the solution Fi is expanded into a B-series, 5(0, xq; h), then a)(0) = 3)„(0) + 
<i>imi&) = 1 implies that (19) can be written as 



5(0), Xo; h) = B((l)„, xq; h) + B(<D,„, B(<D, xq; h); h) 
= Bi^ex, xo; h) + £(0 o (!.,„, xo; h). 



(20) 



Comparing term by term gives us the following result: 
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Theorem 4. The numerical solution Y\ given by{\9) can be written as a B-series 



with <I> recursively defined by 



O(0) = 1, (21a) 
0(T) = 0,,(T) + (O o (D,,„)(T). (21b) 



Since fl>,m(0) = the expression (21b) can be written as 

3>(T) = a)„(T) + a),„(T) + !??(T) 

with 



^(T)= 2 r(T,i9,w) a),>„(,?)]~[cD(,5)l. 



In Table 1 the H functions are listed for aU trees with up to four nodes. 

To decide the weak order we will also need the B-series of the function /, evaluated at the exact and the numerical 
solution. From Theorems 2 and 4 and Lemma 1 we obtain 



f{X{to + h))=Yj m ■ H'Mih) ■ G(m)(xo), 

U€Uf 

f{Yi) = 2 m ■ ipcbium ■ g(m)(xo), 



ueUf 

with 

K 

7=1 

and 

K 

Mm f) = h ^<t(M = [Ti , . . . , rj /)(/l) = W 0(Ty)(/l). 

7=1 

With all the B-series in place, we can now present the order conditions for the weak and strong convergence, for both 
the Ito and the Stratonovich case. As usual we assume that method (1) is constructed such that Eil/^(u){h) - 0(hP^"^) 
"iu € Uf and <1)(t)(/z) = OQi''^'^^) Vr e T, respectively, where especially in the latter expression the C)(-)-notation 
refers to the L^(Q)-norm and h ^ Q. These conditions are fulfilled if for r e T and A: € N = {0, 1, . . . } it holds that 
i^exijyf'' = OOv'-'^P^''^) and (fl>,m(T))^' = Oih^'^P^'^''). Then we have weak consistency of order p if and only if 

EiAo(m)(/i) = Ei/^^(m)(/i) + (9(/i''+') Vm€ t// with p(m) <p+ i (22) 

((22) slightly weakens conditions given in [15]), and mean square global order p if [4] 

^iT)ih) = (piT)ih) + OihP-'h Vre r with p(t) <p, (23) 
E<D(t)(/j) = E(piT)ih) + OihP-'^) Vre T with p(r) </7+ ^ (24) 

and all elementary diff'erentials F(t) fulfill a linear growth condition. Instead of the last requirement it is also enough 
to claim that there exists a constant C such that H^^Cv)!! < C Vy e M", j = 0,...,M, and all necessary partial 
derivatives exist [3]. 
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Table 1: The correction terms H together with the weight functions tp(T) of the correct solution for all trees with up to four nodes. 

T ^(r) (p^ 

Wi 

<l>,>„(.,)0(.,)0(.,) + 0,„, <[)(•«;) + <S>i„ Q . J WjWk ★ dWi 

1^ 3>,m(.,)<D Jl'j + ^im <!>(.,) / / * dWj ★ JW,- 

.A^^^^/ <I>,m(.,)<l>(.;)O(.,)4)(.0 + (!),„. <D(.,)<D(.0 + (D,„. (D(.,)<D(.0 + (D,„. <D(.,)(D(.,) J VK VK ★ JW, 



4. Numerical example 

As an example, we want to construct a parameterized family of strong order 1.5 Taylor methods, applicable to 
Ito-SDEs with one-dimensional noise. As ansatz, we choose 

F„+i =B(a)„, F„; h) + ci7(i)^i,„+i + C2%o,«+i + (^3^(1,1) + C4%'i,„+igi,n+i 

+ C5y^o,«+1^0,n+l + C6 — ^o.n+n^l.n+l'^l.n+l-'' 

where gi„+\ = gK^n+i). ci, ci, C3, C4, C5, e R, 3>eA:(T) = for p(t) > 2, and 3>e;c(T) is required to be deterministic for 
P(t) = 2. 

Together with the order conditions (23) and (24), these requirements are sufficient to uniquely determine the 
remaining unknowns ^ex{T) for p(t) < 2, and we obtain finally the following family of methods, using the abbreviation 

gi = gi(Y„): 

y„+l =y„ + Ci/(i)^i,„+i + C2%0,n+1 + (C3/(l,l) + C4/j)^'i,„+i^l,n+l + C5 y ^o_„+i^O,n+l 

„ 

+ ((-Cl -CA)h + {\ -2ci -C3)I(ij-,)g\gi 

+ ((1 - Ci)/(0,1) - Ci/(i,0))g'ig() + (-<:'2/((),l) + (1 - C2)/(l,0))^o^l 

^Ao,l) - I^Ci + C3 + W(l) + (1 - 3ci - 3c3)/(i,i,i)jg'/(gi,gi) 

- ((Ci + C3 + C4)W(i) + (3C1 + 3C3 - l)/(l,l,l))^'i^'i^l + (1 - 2C2 - €5)^8080 

- (Ci + C4)/l^g'i'(g0,gl) - Clyg'i^o^l - (-C1 + C4)/lVl^'l^0 

+ i(l - 2C2 - C6)h^8oi8u8i) - c^^8\8'x8\8\ - + C3 + C4)h'^8[8ii8u8i) 

3 2 „ , 1 2 ,„ 

-(Cl + -C3 +C4)/l gi{gig,gl)- ^(Ci +C3 +C4)/! {gl,gl,gl). 

In the following, we analyze numerically the order of convergence of this stochastic Taylor methods for the special 
parameter choice ci = C2 = C3 = C4 = C5 = cg = | by applying it to two examples. In both cases, the solution is 
approximated with step sizes 2"^, . . . , 2 and the sample average of M = 2000 independent simulated realizations 
of the absolute error is calculated in order to estimate the expectation. 
Our first example is the non-linear SDE [11] 

Kit) = (iX(5) + yjXis)^ + 1) ds + yjXis^ + ldWis) (25) 

on the time interval / = [0, IJ with the solution X{t) - sinh(f -1- W(t)). 

As second example, we apply our Taylor method to the non-hnear vector-SDE 



w\ ^ r I ix.is) H- vwrpFTTj ^ r /cos(x.(,))\ 

again on the time interval / = [0, 1]. As here we do not know the exact solution, to approximate it we use the 
semi-impUcit order 1.5 method due to Kloeden and Platen [11] with a step size ten times smaller than the actual step 
size. 

For both examples, the results at time ? = 1 are presented in Figure 1, where the order of convergence corresponds 
to the slope of the regression Une. As expected, we observe also numerically strong order 1.5. 
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(a) Error when applied to (25) 



(b) Error when applied to (26) 



Figure 1: Error of example Taylor method vs. step size 



5. Conclusion 

In this article we have now given a new representation formula for (stochastic) B-series evaluated in a B-series, 
which in contrast to the one known for the deterministic case is not based on ordered trees. We have used this formula 
for giving for the first time the order conditions of Taylor methods in terms of rooted trees. We expect these to 
be useful for the construction of new methods e. g. with good stability properties. Finally, we applied these order 
conditions to derive in a simple manner a family of order 1 .5 Taylor methods applicable to Ito SDEs. In a forthcoming 
paper, the theory developed in this article is applied to analyze quite generally the order of iterated Taylor methods 
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